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摘 要 : 针对 中 亚 地 区 的 强生 态 脆 能 性 


高 敏感 性 特征 ,有 必要 开展 广 域 .长 期 的 植被 覆 善 监测 以 匹 


配 "绿色 丝绸 之 路 "的 可 持续 发 展 目标 。 鉴 于 此 ,联合 Landsat 5 ft Landsat 8 卫星 数据 集 , 利 用 Google 


Earth Engine( GEE ) 地 理 空 间 数据 云 计 算 平 台 
度 。 结 果 表 明 :(1) 中 亚 地 区 植被 覆盖 
1993—2018 年间 多 


数 区 域 植被 覆盖 趋势 较为 稳定 ,哈萨克 斯 坦 丘 


,估算 了 中 亚 地 区 1993 一 2018 年 间 共 12 期 的 植被 覆盖 
总 体 水 平 较 低 , 但 也 具有 较为 显著 的 空 


间 异 质 性 ,(2) 中 亚 地 区 
陵 、 费 尔 干 纳 盆地 等 区 域 植被 覆盖 


度 呈 增加 趋势 ,乌拉 尔 河流 域 和 锡 尔 河流 域 等 区 域 植被 覆盖 趋势 为 负 。(3) 植 被 覆盖 度 时 序 特 征 上 ， 
中 亚 地 区 1993 一 2018 年 间 总 体 植被 覆盖 度 累积 增加 3%6, 其 中 吉尔 吉 斯 斯 坦 和 塔吉克 斯 坦 植被 覆盖 


分 别 增加 3.96% 和 5.86%6。(4) 裸 土 区 呈 退 缩 趋势 ,面积 上 
覆盖 区 和 高 植被 覆盖 区 范围 在 呈现 出 的 振荡 式 增 加 。 
能 对 中 亚 地 区 生态 评估 和 演 替 分 析 提 供 技 术 支持 和 定量 


区 进行 区 域 尺 度 的 植被 覆盖 动态 监测 ， 
数据 。 
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中 亚 地 区 地 处 亚 欧 大 陆 腹 地 ,是 典型 的 温 市 沙 
漠 .草原 干旱 区 ,其 生态 环境 承载 力 脆弱 .水 资源 紧 
缺 '，。 近 年 来 ,由 于 人 类 活动 和 全 球 气 候 变 化 对 
中 亚 地 区 生态 环境 的 持续 扰动 ,造成 了 该 区 域 土地 
苑 漠 化 .草原 退化 、 水 资源 过 量 开 发 等 系列 生态 环 
境 问题 ““。 同 时 由 于 特殊 地 缘 条 件 , 中 亚 地 区 被 
BRA -世界 历史 的 地 理 概 组 。 ,也 是 丝绸 之 路 经 济 囊 
的 极为 重要 市 点 之 一 ”。 随 着 “一 带 一 路 ”国家 战略 
的 纵深 推进 ,中 亚 地 区 的 脆弱 生态 环境 将 面临 新 的 
考验 。 为 推行 生态 友好 型 的 发 展 模式 ,从 而 确保 
“绿色 丝绸 之 路 ”可 持续 发 展 , 吸 需 开展 区 域 尺 度 的 
长 时 序 生 态 监 测 和 系统 性 人 研究。 

植被 作为 生态 系统 中 基础 生产 者 ,其 覆盖 状况 
直接 反应 了 生态 系统 的 结构 和 功能 ”。 尤 其 是 以 草 
原生 态 系统 为 主 的 中 亚 地 区 ,植被 覆盖 是 最 为 关键 
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总 计 减 少 25.9x104km2?, 低 植被 覆盖 区 、 中 植被 
研究 结合 遥感 大 数据 和 地 理 云 计算 对 中 亚 地 


植被 覆盖 监测 ; 中 亚 地 区 ; CEE; 多 时 相 Landsat 影 像 ， 归 一 化 植被 指数 
1000 - 6060(2020)04 - 1023 - 10(1023 ~ 1032) 


的 生态 参数 之 一 。 与 传统 的 实地 测 样 、 分 析 的 植被 
监测 方法 相 比 ,遥感 影像 的 植被 监测 手段 的 突出 优 
SET a BI ESSET s aU LCS ES] UNUS S — 
化 植被 覆盖 指数 (Normalized Difference Vegetation 
Index，VDT7) 利 用 植被 对 近 红 外 波段 和 红 光 波段 的 
反射 率 差异 ,能 够 有 效 地 提取 植被 覆盖 信息 , 广 
泛 地 应 用 于 森林 监测 、 草 原 监测 、 农 作物 佑 产 \ 干 
早 监 测 等 诸多 领域 ""“。 因 此 ,基于 NDVI 的 植 
锌 者 盖 估 算 能 够 有 效 测度 区 域 尺 度 的 植被 窗 六 
率 ””。 目 前 ,针对 中 亚 区 域 的 植被 覆盖 遥感 监测 
主要 有 两 方面 :(1) 利 用 中 低 分 辨认 MODerate-reso- 
lution Imaging Spectroradiometer (MODIS) 数据 的 广 
域 植被 覆盖 度 佑 算 及 演变 分 析 研 究 ””。(2) 利 用 
陆地 卫星 (Landsat) 影 像 对 中 亚 局 部 区 域 进行 植被 
T EST. 5 MODIS Wha FA EE , Landsat 卫星 
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的 优势 体现 在 较 高 的 空间 分 辨 率 (30 m) 和 更 长 的 
时 间 路 度数 据 积累 ,能够 为 广 域 .长 时 间 序 列 植被 
监测 提供 更 好 的 数据 文 撑 ”。 但 广 域 、 较 高 分 辩 
x .长 时 间 序 列 的 Landsat 数 据 处 理 、 分 析 受 限于 算 
力 ,目前 还 相当 缺乏 相关 数据 产品 和 人 研究 。GCoogle 
Earth Engine (GEE ) 是 针对 于 地 理 空 间 数据 处 理 、 分 
Tr 可视化 的 云 计算 平台 ,能够 有 效 地 提高 遥感 大 
数据 处 理 的 效率 55 ,是 地 理 大 数据 计算 与 分 析 的 强 
大 工具 。 

鉴于 此 ,顾及 中 亚 区 域 尺度 下 的 数据 空间 覆盖 
率 和 植被 季节 周期 特征 ,利用 GEE 平 台 调 用 和 核查 
了 所 有 可 用 的 Landsat 系 列 卫 星 的 TM 数据 和 OLI 数 
HE , f 1993 4E. . 1994 4E. . 1998 4E. , 2008 Œ , 2009 4E. , 
2010 4F .2013 4F 2014 4F , 2015 4F 2016 4F ,2017 4F 
和 2018 年 等 确定 为 研究 年 份 ,并 估算 了 这 些 年 份 中 
亚 地 区 的 植被 覆盖 度 。 然 后 以 空间 分 布 和 时 序 特 
征 角 度 ,定量 揭示 了 1993 一 2018 年 中 亚 地 区 植被 覆 
盖 度 的 时 空 演 变 规律 。 本 文 能 为 大 尺度 区 域 Land- 
sat 植 被 覆盖 制图 及 其 环境 啊 应 分 析 提 供 可 行 的 技 
术 流 程 和 方法 借鉴 ,同时 长 时 间 序 列 的 植被 履 盖 度 
估算 结果 对 “一 带 一 路 "战略 可 持续 、 深 度 推进 也 具 
有 重要 意义 |。 


1 研究 区 概况 与 数据 概况 


1.1 研究 区 概况 
中 亚 地 区 包括 了 哈萨克 斯 坦 .吉尔吉斯 斯 坦 、 


50° E 55° E 60° E 65° E 


塔吉克 斯 坦 . 乌 效 别克 斯 坦 和 土库曼 斯 坦 等 区 域 。 
该 区 域 的 地 理 位 置 处 于 46°16' E ~ 87°38’ E, 35°11 
N ~55°24' N, 行 政 上 介 于 中 国 、 俄 罗斯 、 伊 明 、 巴 基 
斯 坦 和 阿富汗 等 国之 间 ,面积 共 约 405x10' km"? , 
中 亚 地 区 地 形 呈 现 东南 高 (由 米尔 高 原 . 喀 喇 昆仑 
山脉 、 天 出 山脉 、 阿 尔 泰 山脉 ) 西北 低 ( 图 兰 低 
地 )。 由 于 山脉 和 高 原 阻截 了 来 自 印 度 洋 的 上 暧 湿 气 
流 , 同 时 又 具有 充足 的 日 照 条 件 ,中 亚 地 区 的 蒸发 
量 远 高 于 降水 量 ,从 而 成 为 最 为 典型 的 干旱 区 之 
一 。 区 域内 河流 走向 受到 地 形 影 响 为 西北 走 回 , 15] 
为 无 出 海口 的 内 陆 河 流 。 在 温 市 大 陆 性 气候 的 影 
啊 下 ,中 亚 干旱 区 的 土地 和 窗 被 类 型 主要 是 草地 和 况 
漠 ”。 在 地 缘 条 件 方面 ,中 亚 地 区 是 丝绸 之 路 的 
“ 通 ”" 中 之 重 , 能 源 资源 非常 丰富 ,与 中 国 的 经 济 结 
构 具 有 明显 的 互补 性 ,具有 重要 的 战略 地 位 。 
1.2 数据 概况 

本 人 研究 采用 的 遥感 数据 包括 了 美国 National 
Aeronautics and Space Administration (NASA ) 的 陆地 
卫星 计划 中 Landsat 5 专题 测绘 仪 (Thematic Map- 
per,，TM) 数 据 和 Landsat 8 陆地 成 像 仪 (Operational 
Land Imager, OLD 数据 利用 大 气 辐射 校正 后 的 地 表 
反射 率 产 品 (Land Surface Reflectance, LSR); £8 ii 
中 亚 地 区 的 Landsat 影 像 共 涉及 276 个 条 带 号 /行列 
“+= (28 World Wide Reference System 2, WRS2), 如 
表 1 所 示 。 为 确保 影像 能 够 覆盖 整个 中 亚 地 区 , 采 
用 了 植被 生长 季节 (5.6.7、8. 9 月 ) 所 有 可 用 的 
Landsat 影像 。 由 于 中 亚 区 域 经 纬度 跨度 范 围 大 和 


Ce 


0 250 500 | 
mk | 4 


70? E 75° E 80° E 85° E 


图 1 中 亚 地 区 2018 年 Landsat 8 真 彩色 合成 影像 


Fig. 1 Imagery of the Central Asian countries in 2018 from Landsat 8 true color composite 
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云层 覆盖 影响 ,一 幅 全 区 完整 的 影像 需要 对 大 量 数 
据 进 行 处 理 .滤波 和 拼接 。 如 :2018 年 的 中 亚 地 区 
影像 共处 理 了 2 396 幅 0LI 影 像 。 受 限于 Landsat 5 
TM 影像 的 数据 在 中 亚 地 区 的 时 空 覆 盖 率 ,为 满足 
研究 区 域 数 据 空 间 上 履 善 范围 的 完整 性 ,在 时 间 维 度 
上 ,选择 了 1993 年 .1994 年 1998 4 . 2008 ^F. , 2009 
年 .2010 年 .2013 年 .2014 年 .2015 年 .2016 年 .2017 
年 和 2018 年 为 植被 覆盖 估算 的 年 份 。 其 余年 份 的 
数据 均 无 法 满足 空间 范 于 全 覆盖 的 要 求 。 


表 1 中 亚 地 区 范围 内 Landsat 条 带 号 /行列 号 参照 表 
Tab.1 Reference table of path / row number of Landsat 


within five Central Asian countries 


条 带 号 行列 号 Aa 行列 号 
144 25—27 158 21:393 
145 2327 159 2] 33 
146 25=29 160 22 —34 
147 24-31 161 22 — 34 
148 24 — 32 162 22—34 
149 24 — 32,34 163 22 —34 
150 23 一 34 104 24 — 33 
151 23=33 165 24 — 32 
152 22 —34 166 24 — 30 
133 22 —34 167 24 — 29 
154 22 —34 168 24 — 28 
155 22 — 34 169 24 — 28 
156 2] 35 170 24 - 27 
157 21-36 171 25-26 


2 研究 方法 


2.1 数据 预 处 理 

不 同 于 单 景 影像 预 处 理 流程 ,GEE 以 数据 集 概 
念 切入 建立 数据 集 预 筛选 .处 理工 作 流 。 由 于 直接 
调用 了 GEE 中 的 LSR 产 品 (LEDAPS 校 正 )” ,而 不 
再 需要 辐射 校正 处 理 ,技术 路 线 如 图 2 所 示 。 数 据 
筛选 的 条 件 :(1) 拍摄 覆盖 包含 中 亚 地 区 的 所 有 影 
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像 .(2) 观测 时 间 为 植被 生长 季 市 (5~9 月 份 ) 期 间 
的 所 有 影像 。 数 据 预 处 理 主要 包括 :(1) 利用 CF- 
MASK 算 法 进行 云 BRE EE" (2) 采用 时 序 上 的 
中 值 滤 波 方 法 使 每 年 的 数据 集 融 合 为 一 景 范围 全 
和 覆盖、 质量 均衡 平稳 的 无 云 影像 .(3) 影像 裁剪 HE 
接 。(4) 逐一 核查 年 合成 影像 ,剔除 由 于 Landsat 数 
据 集 本 身 空间 覆盖 率 不 足 造成 的 空间 覆盖 率 低 的 
年 合成 影像 。 

中 亚 地 区 Landsat 影 像 缺 失 现 象 较为 明显 ,因此 
通过 增加 数据 采集 时 间 ,来 提高 大 尺度 下 的 影像 空 
间 有 覆盖 率 。 同 时 ,通过 时 序 上 中 值 滤 波 处 理 可 以 有 
效 地 减 小 不 同 条 带 号 由 于 观测 期 不 同 造成 的 影像 
色差 ,增强 了 影像 质量 的 稳定 性 。 最 后 ,不同 年份 
的 数据 集 进 行 中 值 滤波 处 理 后 ,提高 了 影像 在 时 序 

上 的 鲁 棒 性 和 可 对 比 性 。 
2.2 植被 覆盖 估算 

归 一 化 植被 指数 是 遥感 影像 中 近 红 外 波段 的 反 
射 值 和 红 光 波段 的 反射 值 之 差 与 其 之 和 的 比值 ( 式 
1) ,是 反映 植被 覆盖 的 重要 参数 之 一 。Landsat 5 TM 
影像 的 3、4 波 段 分 别 对 应 了 红 光 波段 和 近 红 外 波 
段 ,而 Landsat 8 OLI 影 像 的 红 光 波段 和 近 红 外 波段 
分 别 式 4.5 波 段 。 

NDVI= ELDER (1) 


像 元 二 分 模型 假设 地 表单 元 是 由 植被 覆盖 与 
非 植被 覆盖 构成 ,是 混合 像 元 分 解 模型 中 最 简单 的 
模型 。 但 其 形式 简洁 而 且 具 有 一 定 物 理 意义 ,因此 
广泛 应 用 于 植被 覆盖 度 估算 研究 中 。 针 对 中 亚 
地 区 的 植被 履 盖 估算 ,已 有 学 者 分 析 了 像 元 二 分 模 
型 的 适用 性 和 可 靠 性 后 。 鉴 于 此 ,采用 像 元 二 分 模 
型 ( 式 2) 估 算 中 亚 地 区 的 植被 覆盖 度 : 

NDVI - NDVI „ 

NDVI,- NDVI 
式 中 :FVC 是 植被 覆盖 度 ;NDVI 是 归 一 化 植被 指数 ; 
NDVLw 是 裸 土 的 NDVI 什 ;NDVI 是 全 植被 像 元 的 


FVC- (2) 


图 2 基于 CEE 平 台 LSR 数 据 筛选 及 预 处 理工 作 流 
Fig.2 The workflow of data filtering and preprocessing of LSR based on CEE platform 
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NDVI. PEHE p BS FE RC s Pec en HJ EX Sak Op f TE 
天 山 山 脉 .阿尔 泰山 山脉 等 区 域 , 此 类 区 域 得 益 于 
雪山 融化 补给 地 表 水 ,覆盖 着 西伯 利 亚 落 叶 松 和 和 雪 
岭 云 杉 等 针 叶 林 ,植被 履 盖 最 大 的 接近 于 1, 其 NDVI 
多 年 均值 约 为 0.9。 卡 拉 库 姆 沙漠 、 克 孜 勒 库 姆 沙漠 
等 沙漠 区 域 无 植被 履 盖 ,此 类 区 域 的 NDVI<0.1。 
此 ,本 文中 NDVL 设 为 0.9,NDVL 设 为 0.1。 

采用 一 元 线性 回归 分 析 ,分 析 植 被 覆盖 度 序列 
数据 趋势 。 趋 势 分 析 对 每 个 像 元 的 植被 覆盖 度 进 
行 时序 建 模 ,通过 其 变化 特征 反应 一 定时 间 序 列 的 
植被 覆盖 格局 演变 。 其 计算 公式 如 下 : 


n YY,x FVC, - VY, V. FVC, 
_ ii T Get 


(3) 


式 中 :0 是 植被 覆盖 度 FVC; 回 归 方 程 的 斜率 ;n 是 
植被 监测 年 数 ;是 序号 i 对 应 的 年 份 。 若 04w 值 为 
正 , 则 表明 植被 覆 次 度 增加 ,反之 表明 植被 履 兰 
退化 。 


3 结果 与 分 析 


3.1 植被 覆盖 度 估算 结果 
中 亚 地 区 1993 一 2018 年 间 植 被 覆盖 测度 结 
如 图 3 所 示 , 其 中 水 域 (湖泊 、 河 流 以 及 积 雪 冰川 ) 是 
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由 NDVI 负 值 区 域 掩 膜 而 来 。 参 考 等 间距 原则 并 结 
合 Google Earth 高 分 影像 采集 的 样本 ,将 中 亚 地 区 植 
被 覆盖 测度 结果 划分 为 以 下 4 个 等 级 :(1) 裸 土 , 即 
无 植被 覆盖 区 域 (FYC<0) 。(2) 低 植被 覆盖 区 域 
(0<FYC< 和 0.3)。(3) 中 植被 覆盖 区 域 (0.3 < FVC < 
0.6)。(4) 高 植被 覆盖 区 域 (0.6 < FVC 1). 

在 总 体 上 ,中 亚 地 区 整体 的 植被 覆盖 的 空间 
格局 表现 为 植被 覆盖 度 较 低 (2018 年 区 域 植被 覆 
盖 度 为 0.18) ,反应 了 中 亚 地 区 大 部 分 区 域 是 受 干 
燥 气 候 条 件 影 响 而 形成 的 生态 脆弱 干旱 区 。 但 是 
植被 覆盖 的 空间 格局 也 存在 异 质 性 ,主要 有 以 下 
三 个 典型 区 域 :(1) 图 兰 平原 的 植被 呈现 出 大 范围 
的 低 覆 盖 。 尤 其 是 该 区 的 沙漠 区 域 ( 占 平原 面积 
一 半 ) 植 被 覆盖 极为 稀 玻 趋 近 于 无 植被 覆盖 ,其 中 
主要 包括 卡拉 库 姆 沙漠 (土库曼 斯 坦 境内 ) , vd 
勒 库 姆 沙漠 (地 路 乌 效 别克 斯 坦 、 哈 萨克斯 坦 )。 
(2) 中 亚 地 区 境内 由 东南 各 东 北方 回 分 布 有 帕 米 
尔 高 原 、 喀 拉 昆 仑 山脉 .昆仑 山脉 、 天 山 山 脉 和 阿 
尔 泰山 脉 ,由 于 夏季 高 原 冰川 融化 ,水 资源 补给 较 
足 , 使 这 些 区 域 的 植被 覆盖 情况 较 好 ,在 图 3 中 呈 
现 出 浅 绿 色 和 深 绿 色 。(3) 中 亚 区 域 北 部 的 哈萨克 
斯 坦 丘 陵 区 域 形 成 迎风 坡 ,阻挡 来 自 大 西洋 的 暖 
湿 气 流 从 而 形成 降雨 ,因此 该 区 域 植被 覆盖 情况 
也 相对 较 好 。 
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图 3 中 亚 地 区 1993 一 2018 年 植被 覆盖 度 空 间 分 布 


Fig. 3 Spatial distribution of vegetation coverage over the Central Asian countries between 1993 and 2018 


3.2 植被 覆盖 空间 格局 演变 分 析 
中 亚 地 区 1993 一 2018 年 植被 覆盖 度 趋势 变化 
空间 格局 如 图 4 所 示 ( 图 中 空白 区 域 是 由 历年 Land- 
sat 影像 缺失 区 累积 而 成 )。 从 植被 覆盖 变化 趋势 空 
间 格 局 上 看 ,中 亚 地 区 1993 一 2018 年 期 间 大 部 分 区 
域 植 被 覆盖 度 保持 在 较 稳定 状态 。 图 4 中 ,植被 柳 
盖 度 呈 增 加 的 趋势 的 区 域 主要 有 :(1) 天 山 山 脉 由 
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漠 化 等 植被 覆盖 趋势 下 降 的 区 域 主要 在 马 拉 尔 河 
流域 (C 区 域 ) 以 及 位 于 到 和 孜 勒 库 姆 沙漠 以 北 锡 尔 河 
流域 (D 区 域 )。 同 时 ,可 以 发 现 原 世界 第 四 大 淹 的 
威海 区 域 (A 区 域 ) 的 植被 覆盖 测度 值 变化 趋势 较为 
明显 ,但 需要 说 明 的 是 该 现象 并 不 是 由 植被 覆盖 度 
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增加 造成 , 而 是 土地 覆 被 类 型 由 水 体 转移 为 裸 土 导 
致 。 戌 海流 域 由 于 过 量 的 工农 业 用 水 导致 其 水 位 
下 降 面积 萎缩 ,从 而 裸露 出 湖 床 沉积 层 。 

采用 下 检验 对 中 亚 地 区 1993 一 2018 年 的 植被 
履 盖 趋势 进行 显著 性 检验 ,检验 后 的 结果 被 分 为 五 
类 。 其 中 ,变化 不 显著 (P>0.05) 区 域 面积 ， 
348.56 x 10* km’, 占 比 86.03% ;显著 减少 (0,,,<0， 
0.01 < Px 0.05) E fi t E LIP (Ope < 0, P<O.01) K 
Ty rfi FALE 4.85 x 10* km’, E EE 1.199% ;显著 增加 
(bawe > 0, 0.01 < P 0.05) 与 极 显著 增加 (be> 0, 
P<0.01) HRF W 18.88x 10* km’, H HE 5.07%., 
EHE X EA SHEL BZ e 9425] (8.0.0029 ,表明 了 中 亚 
地 区 总 体 较 为 稳定 的 植被 覆盖 变化 趋势 。 如 表 2 所 
示 , 包 含 卡 拉 库 姆 沙漠 的 土库曼 斯 坦 植被 覆盖 趋势 
测度 值 最 小 ;地 处 天 山 山 脉 西 段 的 吉尔 吉 斯 斯 坦 的 
植被 覆盖 多 样 性 和 复杂 性 较 强 ,植被 覆盖 趋势 增 
加 ,但 由 于 该 区 域 植 被 空间 分 布 异 质 性 强 ,造成 该 
区 域 植 被 覆盖 趋势 的 标准 差 较 大 ;土库曼 斯 坦 境内 
土地 利用 / 复 盖 类 型 较为 单一 ,主要 由 沙漠 和 裸 土地 
表 构 成 ,该 区 域 空 间 同 质 性 强 , 造 成 植被 覆盖 变化 


35?N 
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1998 一 2018 年 ,中 亚 地 区 植被 覆盖 度 成 振荡 式 增 加 
趋势 ,期 间 植 被 履 盖 度 累 积 增加 10%。 

中 亚 地 区 1993 年 .1994 年 .1998 年 和 2008 年 植 
被 覆盖 度 空间 分 布 较为 稳定 ,其 标准 差 均 低 于 0.3。 
但 此 之 后 ,植被 覆盖 度 在 空间 上 趋 于 分 化 ,吉尔 吉 
斯 斯 坦 和 塔吉克 斯 坦 的 空间 异 质 性 尤为 显著 。 上 


表 2 中 亚 地 区 1993 一 2018 年 植被 覆盖 趋势 统计 表 
Tab.2 Statistic table of vegetation coverage trend over 
the Central Asian countries between 1993 and 2018 


国家 Slope 均值 /103 Slope 标准 差 / 107 
吉尔 吉 斯 斯 坦 3.50 16.02 
II D vC SH 3.23 4.66 
塔吉克 斯 坦 3.12 8.10 
土库曼 斯 坦 0.99 3.25 
乌兹别克 斯 坦 1.97 5.87 


3.3 不 同 国家 植被 覆盖 度 时 序 演变 分 析 
1993—2018 年 中 亚 地 区 植被 覆盖 度 均值 变化 
曲线 如 图 5(a) 所 示 ,中 亚 地 区 植被 覆盖 度 由 高 至 低 
依次 是 吉尔 吉 斯 斯 坦 .哈萨克 斯坦 .塔吉克 斯 坦 S 
兹 别克 斯 坦 和 土库曼 斯 坦 。 总 体 上 ,中 亚 地 区 
1993—2018 年 间 ,植被 覆盖 度 累积 增加 3%。 其 中 ， 
植被 禾 盖 累积 增加 率 最 高 的 是 塔吉克 斯 坦 ,达到 
5.86%; 而 土库曼 斯 坦 由 于 境内 几 无 植被 覆盖 ,25 a 
间 植 被 覆盖 度 基 本 不 变 。1993 一 1998 年 期 间 ,中 亚 
地 区 植被 覆盖 总 体 呈 显著 减 小 趋势 ,达到 — 7% ; 而 
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ik: A 威海 区 域 ; B 哈萨克 丘陵 区 域 ; C 乌拉 尔 河流 域 ; D 锡 尔 河流 域 
图 4 中 亚 地 区 1993 一 2018 年 植被 覆 兰 度 变 化 趋势 空间 格局 

Fig.4 Spatial pattern of vegetation coverage trend over the Central Asian 


countries between 1993 and 2018 


述 两 国境 内 坐落 有 天 山 山 脉 和 由 米尔 高 原 ,区 域内 
的 冰川 雪山 对 全 球 气候 变化 具有 较 强 的 敏感 性 ,从 

e Je i^ Ce D ,使 植被 覆 兰 度 及 其 空间 分 布 呈 
现 出 较 强 的 时 序 波动 特征 。 研 究 结 果 与 利用 中 分 
辩 率 MODIS 数 据 估算 中 亚 地 区 草地 覆盖 度 时 空 特 
征 保持 较 高 的 一 致 性 ”。 低 分 辨识 GIMMS 数 据 的 
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图 5 ”中 亚 地 区 1993 一 2018 年 植被 覆盖 度 时 序 演变 特征 


Fig.5 Temporal features of vegetation coverage over the Central Asian countries between 1993 and 2018 


估算 结果 表明 ,中 亚 地 区 1982—2012 年 期 间 NDVI 
呈现 出 微弱 增长 趋势 " ,也 吻合 本 文 结论 。 
3.4 不 同 植被 覆盖 区 面积 时 序 变化 分 析 

通过 空间 统计 各 期 影像 1993 一 2018 年 水 域 . 裸 
土 、 低 植被 覆盖 .中 植被 覆盖 和 高 植被 覆盖 面积 ,其 
占 比 结果 如 图 6 所 示 。 总 体 上 ,1993 一 2018 年 监测 
期 间 ,无 植被 覆盖 裸 土 的 覆盖 范围 呈 退 缩 趋 势 , 面 
TUTUP 25.9x10' km:。 低 植被 覆盖 区 总 体 上 呈 
现 出 扩张 趋势 ,1993 一 2018 年 其 面积 总 计 增 加 
19.6x10' km。 中 植被 覆盖 区 和 高 植被 覆盖 区 面积 
在 1993 一 2018 年 间 振 沪 式 增加 , 累积 增加 量 分 别 为 
9.5x10* km? All 7.1x 10° km。 水 域 主要 由 于 成海 收 
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图 6 中 亚 地 区 1993 一 2018 年 不 同 区 域 植被 覆盖 面积 占 比 
变化 时 序 特 征 
Fig. 6 Temporal features of area proportion of different regions 


within the Central Asian countries between 1993 and 2018 
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HÄ , 1993—2018 年 间 面 积 总 计 减 少 3.3x10' km’, H 
中 1998 - 2008 年 减少 量 最 为 剧烈 ,减少 面积 达到 
4.3x10* km , 占 累积 损失 量 的 82% ,表明 了 此 期 间 中 
亚 地 区 干燥 气候 对 地 表 环 境 造 成 较 强 的 扰动 。 同 
时 ,1998 年 一 2008 年 间 低 植被 覆盖 面积 总 共 减 少 
38.5x10* km ,而 裸 土 面积 增加 了 28.1km ,也 反映 了 
该 期 间 区 域 生态 环境 的 退化 较为 明显 。 


4 结论 


2013 年 “一 带 一 路 ”的 国家 战略 ( 共 建 “丝绸 之 
路 经 济 带 ” 和 “21 世纪 海上 丝绸 之 路 ”合作 倡议 BS 
提出 使 中 亚 地 区 的 地 理 区 域 优势 极 具 突出 。 中 亚 
地 区 作为 “一 带 一 路 ”的 重要 节点 和 通道 ,决定 了 中 
国 与 其 的 联系 将 更 加 紧密 从 而 形成 “命运 共同 
体 ”"。 倡 导 既 要 “金山 ”“ 银 山 ”, 又 要 绿 水 青山 的 发 
展 模式 对 区 域 经 济 结构 和 生态 环境 的 协调 匹配 度 
提出 了 更 高 .更 科学 的 要 求 。 在 此 ,本 研究 采用 
Landsat 长 时 序数 据 , 利 用 GEE 云 计算 平台 ,实现 了 
中 亚 地 区 1993—2018 年 区 域 大 尺度 、 较 高 分 辩 率 
(30 m) 的 植被 覆盖 度 佑 算 ,并 以 空间 格局 和 时 序 特 
征 角度 切入 具体 分 析 了 中 亚 地 区 植被 覆盖 的 时 空 
演变 规律 。 结 果 表 明 : 

(1) 在 空间 格局 上 ,中 亚 区 域 植被 覆盖 具有 较 
强 的 空间 异 质 性 ,其 典型 植被 覆盖 区 域 包含 :由 低 
植被 覆盖 荒漠 和 无 植被 覆盖 沙漠 构成 的 图 兰 低 地 ， 
较 高 植被 覆盖 的 哈萨克 斯 坦 丘 陵 以 及 高 植被 覆盖 
的 东部 天 山 、 阿 尔 泰山 区 域 。 
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(2) 植 被 覆盖 趋势 空间 特征 上 ,中 亚 地 区 
1993 一 2018 年 间 多 数 区 域 植 被 覆盖 趋势 较为 稳定 ， 
而 植被 覆盖 趋势 增加 主要 分 布 在 东部 费 尔 干 那 盆 
地 哈萨克 斯 坦 中 部 丘陵 区 域 ,趋势 减少 的 区 域 主 
要 有 乌拉 尔 河 流域 和 锡 尔 河流 域 。 

(3) 植被 覆盖 度 时 序 特征 上 ,中 亚 地 区 1993 一 
2018 年 间 总 体 植被 覆盖 度 累积 增加 3%。 吉 和 尔 吉 其 
斯 坦 和 塔吉克 斯 坦 植被 覆盖 呈现 出 较为 显著 的 增 
加 趋势 ,分 别 增加 3.96% 和 5.86%。 土 库 曼 斯 坦 境内 
多 为 荒漠 和 沙漠 , 其 植被 覆盖 常年 维持 不 变 。 

(4) 植被 禾 盖 类 别 范 围 演变 特征 上 , 1993 一 
2018 年 无 植被 覆盖 裸 土 的 覆盖 范围 呈 退 缩 趋势 , 面 
只 总 计 减 少 25.9x10' km。 水 域 主要 由 于 成 海 收缩 ， 
面积 总 计 减 少 $.3x10' km?。 低 植被 覆盖 区 .中 植被 
覆盖 区 和 高 植被 覆盖 区 范围 在 1993 一 2018 年 期 间 
呈现 出 的 振荡 式 扩张 特征 。 

综 上 所 述 ,本 文通 过 利用 Landsat 遥感 数据 , 结 
合 GEE 平 台 , 从 空间 大 尺度 和 较 高 分 辨 率 上 提取 分 
析 中 亚 地 区 的 植被 覆盖 情况 以 及 演变 规律 。 该 研 
究 结果 有 利于 从 整体 上 了 解 中 亚 地 区 的 植被 覆盖 
状况 演变 规律 乃至 湖泊 水 面 变化 规律 ,从 而 为 "一 
带 一 路 "国家 战略 的 精准 、 纵 深 推 进 提供 空间 大 尺 
度 的 生态 环境 参考 和 数据 基础 。 
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Abstract: Central Asian regions with a typical temperate desert as well as grassland environment experience sev- 
eral ecological environmental problems such as water shortages and land degradation. Particularly in recent de- 
cades, the ecological disturbances caused by climate change and human activity have led to more consideration be- 
ing given to environmental protection. As the most essential producer of an ecosystem, vegetation coverage status 
can directly affect an ecosystem's function and structure. Therefore, to align with the sustainable development goal 
of the "Green Silk Road" ,there is a need to conduct vegetation coverage monitoring with wide-range and long-term 
abilities. In view of this, 12 periods of vegetation coverage in the Central Asian regions from 1993 to 2018 were esti- 
mated using the Google Earth Engine geospatial data cloud computing platform from Landsat 5 together with Land- 
sat 8 remote sensing data sets. To start, the research area boundary was used as the selected condition to filter the 
Landsat images. Through data preprocessing, including a cloud mask, filter, and split, we acquired 12 high-quality 
remote sensing images over the Central Asian regions. Then, we calculated the normalized difference vegetation in- 
dex using the NIR and RED bands of the selected images. The dimidiate pixel model was used to estimate the vege- 
tation coverage , followed by a linear regression analysis. The results showed as follows: (1) The overall level of veg- 
etation coverage in Central Asia was low but had significant spatial heterogeneity. (2) The vegetation coverage 
change trend of most regions in Central Asia between 1993 and 2018 was relatively stable. The vegetation coverage 
over the Kazakhstan Hills and Fergana Valley areas showed an increasing trend, while the vegetation coverage in 
the Ural River Basin and Syr Darya River Basin areas showed a negative trend. (3) In terms of the time series char- 
acteristics of the vegetation coverage, the total vegetation coverage increased by 3% in the Central Asian regions be- 
tween 1993 and 2018, during which the vegetation coverage in Kyrgyzstan and Tajikistan increased by 3.96% and 
5.86% , respectively. (4) The bare soil area showed a retreating trend, and its total area decreased by 259 000 km’. 
(5) Lastly, the low vegetation, middle vegetation, and high vegetation cover areas showed an oscillatory increase. 
Overall, this study combined remote sensing data with geographical cloud computing to monitor vegetation coverage 
on a regional scale in Central Asia to provide technical support and quantitative data for ecological assessment and 
succession analysis in Central Asia. 

Key words: vegetation coverage monitoring; the Central Asian Regions; Google Earth Engine; multi-temporal 
Landsat images ; NDVI 


